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This review is written at the time of the twentieth anniversary of the discovery of high temperature super- 
conductors, which, nearly coincides with the important discovery of the superfiuid phases of ultracold trapped 
fermionic atoms. We show how these two subjects have much in common. Both have been addressed from the 
perspective of the BCS-Bose Einstein condensation (BEC) crossover scenario, which is designed to treat short 
coherence length superfluids with transition temperatures which are "high", with respect to the Fermi energy. A 
generalized mean field treatment of BCS-BEC crossover at general temperatures T, based on the BCS-Leggett 
ground state, has met with remarkable success in the fermionic atomic systems. Here we summarize this suc- 
cess in the context of four different cold atom experiments, all of which provide indications, direct or indirect, 
for the existence of a pseudogap. This scenario also provides a physical picture of the pseudogap phase in the 
underdoped cuprates which is a central focus of high T c research. We summarize successful applications of 
BCS-BEC crossover to key experiments in high T c systems including the phase diagram, specific heat, and 
vortex core STM data, along with the Nernst effect, and exciting recent data on the superfiuid density in very 
underdoped samples, 
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in the field of high T c superconductivity 
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and its implications below T c . Either the excitation gap is inti- 
mately connected to the superconducting order reflecting, for 
example, the existence of "pre-formed pairs", or it is extrin- 
sic and associated with a competing ordered state unrelated to 
superconductivity. 

The emphasis of this Review is on the pseudogap state as 
addressed by a particular preformed pair scenario which has 
its genesis in what is now referred to as "BEC-Bose Einstein 
condensation (BEC) crossover theory". Here one contem- 
plates that the attraction (of unspecified origin) which leads 
to superconductivity is stronger than in conventional super- 
conductivity. In this way fermion pairs form before they Bose 
condense, much as in a Bose superfluid. In support of this 
viewpoint for the cuprates are the observations that: (i) the co- 
herence length £ for superconductivity is anomalously short, 
around 10A as compared with 1000A for a typical supercon- 
ductor. Moreover (ii) the transition temperatures are anoma- 
lously high, and (iii) the systems are close to two dimen- 
sional (2D) (where pre-formed pair or fluctuation effects are 
expected to be important). Finally, (iv) the pseudogap has the 
same d-wave symmetry IU6I1 as the superconducting order pa- 
rameter y|,|4j] and there seems to be a smooth evolution of the 
excitation gap from above to below T c . 

To investigate this BCS-BEC crossover scenario we have 
the particular good fortune today of having a new class 
of atomic physics experiments involving ultracold trapped 
fermions which, in the presence of an applied magnetic field, 
have been found to have a continuously tunable attractive in- 
teraction. At high fields the system exhibits BCS-like super- 
fluidity, whereas at low fields one sees BEC-like behavior. 

This Review presents a consolidated study of both the pseu- 
dogap phase of the cuprates and recent developments in ultra- 
cold fermionic superfluids. The emphasis of these cold atom 
experiments is on the so-called unitary or strong scattering 
regime, which is between the BEC and BCS limits, but on 
the fermionic side. The superfluid state in this intermediate 
regime is also referred to in the literature as a "resonant super- 
fluid" HlTlflSll . Here we prefer to describe it as the "pseudogap 
phase", since that is more descriptive of the physics and un- 
derlines the close analogy with high T c systems. Throughout 
this Review we will use these three descriptive phrases inter- 
changeably. 



Fermionic Pseudogaps and Meta-stable Pairs: Two Sides of 
the Same Coin 



BCS-BEC crossover theory is based on the observations of 
Eagles 1 19] and Leggett 1 20] who independently noted that 
the BCS ground state wavefunction 



¥o=n k (u k + «k4<£ k )|0) 



(1) 



had a greater applicability than had been appreciated at the 
time of its original proposal by Bardeen, Cooper and Schrief- 
fer (BCS). As the strength of the attractive pairing interaction 
U (< 0) between fermions is increased, this wavefunction is 
also capable of describing a continuous evolution from BCS 
like behavior to a form of BEC. What is essential is that the 



chemical potential fj, of the fermions be self consistently com- 
puted as U varies. 

The variational parameters u k and u k are usually repre- 
sented by the two more directly accessible parameters A sc (0) 
and /i, which characterize the fermionic system. Here A sc (0) 
is the zero temperature superconducting order parameter. 
These fermionic parameters are uniquely determined in terms 
of U and the fermionic density n. The variationally deter- 
mined self consistency conditions are given by two BCS-like 
equations which we refer to as the "gap" and "number" equa- 
tions respectively. 



A sc (0) 
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2E k 



where 



E k = V(£k - M) 2 + A2 C (0) 



(2) 



(3) 



and e k = h 2 k 2 /2m are the dispersion relations for the Bo- 
goliubov quasiparticles and free fermions, respectively. An 
additional advantage of this formalism is that Bogoliubov de 
Gennes theory, a real space implementation of this ground 
state, can be used to address the effects of inhomogeneity and 
external fields at T = 0. This has been widely used in the 
crossover literature. 

Within this ground state there have been extensive studies 
l2lll of collective modes l22L l23ll and effects of two dimen- 
sionality l22ll . Nozieres and Schmitt-Rink were the first l24ll 
to address non-zero T. We will outline some of their con- 
clusions later in this Review. Randeria and co-workers refor- 
mulated the approach of Nozieres and Schmitt-Rink (NSR) 
and moreover, raised the interesting possibility that crossover 
physics might be relevant to high temperature superconduc- 
tors l22ll . Subsequently other workers have applied this pic- 
ture to the hi gh T r cuprates 12511211 l27ll and ultracold fermions 
1171 Il8l 1281 129II as well as formulated alternative schemes 
I30t Bill for addressing T 7^ 0. Importantly, a number of ex- 
perimentalists, most notably Uemura 02l> have claimed evi- 
dence in support I33ll34t 13511 of the BCS-BEC crossover pic- 
ture for high T c materials. 

Compared to work on the ground state, considerably less 
has been written on crossover effects at non-zero temperature 
based on Eq. Q. Because our understanding has increased 
substantially since the pioneering work of NSR, and because 
they are the most interesting, this review is focused on these 
finite T effects. 

The importance of obtaining a generalization of BCS the- 
ory which addresses the crossover from BCS to BEC ground 
state at general temperatures T < T c cannot be overesti- 
mated. BCS theory as originally postulated can be viewed 
as a paradigm among theories of condensed matter systems; it 
is complete, in many ways generic and model independent, 
and well verified experimentally. The observation that the 
wavefunction of Eq. Q goes beyond strict BCS theory, sug- 
gests that there is a larger mean field theory to be addressed. 
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FIG. 1: Behavior of the T = chemical potential /i in the three 
regimes, /i is essentially pinned at the Fermi temperature Ef in the 
BCS regime, whereas it becomes negative in the BEC regime. The 
PG (pseudogap) case corresponds to non-Fermi liquid based super- 
conductivity in the intermediate regime. 



Equally exciting is the possibility that this mean field the- 
ory can be discovered and simultaneously tested in a very 
controlled fashion using ultracold fermionic atoms fnl [Till . 
Mean field approaches are always approximate. We can as- 
cribe the simplicity and precision of BCS theory to the fact 
that in conventional superconductors the coherence length £ is 
extremely long. As a result, the kind of averaging procedure 
implicit in mean field theory becomes nearly exact. Once £ 
becomes small BCS is not expected to work at the same level 
of precision. Nevertheless even when they are not exact, mean 
field approaches are excellent ways of building up intuition. 
And further progress is not likely to be made without investi- 
gating first the simplest of mean field approaches, associated 
with Eq. Q. 

The effects of BEC-BCS crossover are most directly re- 
flected in the behavior of the fermionic chemical potential p. 
We plot the behavior of p in Fig.^ which indicates the BCS 
and BEC regimes. In the weak coupling regime /i = Ep and 
ordinary BCS theory results. However at sufficiently strong 
coupling, /i begins to decrease, eventually crossing zero and 
then ultimately becoming negative in the BEC regime, with 
increasing \U\. We generally view /i = as a crossing point. 
For positive p the system has a remnant of a Fermi surface, 
and we say that it is "fermionic". For negative p, the Fermi 
surface is gone and the material is "bosonic". 

The new and largely unexplored physics of this problem lies 
in the fact that once outside the BCS regime, but before BEC, 
superconductivity or superfluidity emerge out of a very exotic, 
non-Fermi liquid normal state. Emphasized in Fig. ^ is this 
intermediate (i.e., pseudogap or PG) regime having positive 
p which we associate with non-Fermi liquid based supercon- 
ductivity 1 251 \M& 13711 . Here, the onset of superconductivity 
occurs in the presence of fermion pairs. Unlike their counter- 
parts in the BEC limit, these pairs are not infinitely long lived. 
Their presence is apparent even in the normal state where an 
energy must be applied to create fermionic excitations. This 
energy cost derives from the breaking of the metastable pairs. 
Thus we say that there is a "pseudogap" (PG) at and above 
T c . It will be stressed throughout this Review that gaps in the 
fermionic spectrum and bosonic degrees of freedom are two 
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FIG. 2: Contrasting behavior of the excitation gap A(T) and or- 
der parameter A„ C (T) versus temperature in the pseudogap regime. 
The height of the shaded region reflects the number of noncondensed 
pairs, at each temperature. 



sides of the same coin. A particularly important observation 
to make is that the starting point for crossover physics is based 
on the fermionic degrees of freedom. A non-zero value of the 
excitation gap A is equivalent to the presence of metastable 
or stable fermion pairs. And it is only in this indirect fashion 
that we can probe the presence of these "bosons", within the 
framework of Eq. Q. 

In many ways this crossover theory appears to represent a 
more generic form of superfluidity. Without doing any calcu- 
lations we can anticipate some of the effects of finite tempera- 
ture. Except for very weak coupling, pairs form and condense 
at different temperatures. More generally, in the presence of 
a moderately strong attractive interaction it pays energetically 
to take some advantage and to form pairs (say roughly at tem- 
perature T*) within the normal state. Then, for statistical rea- 
sons these bosonic degrees of freedom ultimately are driven 
to condense at T c < T*, as in BEC. 

Just as there is a distinction between T c and T*, there must 
be a distinction between the superconducting order param- 
eter A sc and the excitation gap A. In Fig. [2] we present a 
schematic plot of these two energy parameters. It may be seen 
that the order parameter vanishes at T c , as in a second order 
phase transition, while the excitation gap turns on smoothly 
below T*. It should also be stressed that there is only one 
gap energy scale in the ground state 1 20] of Eq. Q. Thus 
A sc (0) = A(0). 

In addition to the distinction between A and A sc , another 
important way in which bosonic degrees of freedom are re- 
vealed is indirectly through the temperature dependence of A. 
In the BEC regime where fermionic pairs are pre-formed, A 
is essentially constant for all T < T c (as is By contrast in 
the BCS regime it exhibits the well known temperature depen- 
dence of the superconducting order parameter. This behavior 
is illustrated in Fig. [3] 

Again, without doing any calculations we can make one 
more inference about the nature of crossover physics at finite 
T. The excitations of the system must smoothly evolve from 
fermionic in the BCS regime to bosonic in the BEC regime. 
In the intermediate case, the excitations are a mix of fermions 
and meta-stable pairs. Figure |4] characterizes the excitations 
out of the condensate as well as in the normal phase. This 
schematic figure will play an important role in our thinking 
throughout this Review. 
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FIG. 3: Comparison of temperature dependence of excitation gaps 
in BCS (left) and BEC (right) limits. The gap vanishes at T c for the 
former while it is essentially T-independent for the latter. 
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FIG. 4: The character of the excitations in the BCS-BEC crossover 
both above and below T c . The excitations are primarily fermionic 
Bogoliubov quasiparticles in the BCS limit and bosonic pairs (or 
"Feshbach bosons") in the BEC limit. For atomic Fermi gases, the 
"virtual molecules" in the PG case consist primarily of "Cooper" 
pairs of fermionic atoms. 



Introduction to high T c Superconductivity: Pseudogap 
Effects 



This Review deals with the intersection of two fields and 
two important problems: high temperature superconductors 
and ultracold fermionic atoms in which, through Feshbach 
resonance effects, the attractive interaction may be arbitrarily 
tuned by a magnetic field. Our focus is on the broken symme- 
try phase and how it evolves from the well known ground state 
at T = to T = T c . We begin with a brief overview JHH 
of pseudogap effects in high temperature superconductors. A 
study of concrete data in these systems provides a rather nat- 
ural way of building intuition about non-Fermi liquid based 
superfluidity, and this should, in turn, be useful for the cold 
atom community. 

It has been argued by some jHHHllllES that a BCS- 
BEC crossover-induced pseudogap is the origin of the mys- 
terious normal state gap observed in high temperature super- 
conductors. While this is a highly contentious subject some of 
the arguments in favor of this viewpoint (beyond those listed 
in Section ITXl rest on the following observations: (i) To a 
good approximation the pseudogap onset temperature lEHI^ 
T* fn 2A(0)/4.3 which satisfies the BCS scaling relation, 
(ii) There is widespread evidence for pseudogap effects both 
above 0] as well as (iii) below j4j T c . (iv) In addi- 
tion, it has also been argued that short coherence length su- 
perconductors may quite generally exhibit a distinctive form 
of superconductivity [32] which sets them apart from conven- 
tional superconductors. One might want, then, to concentrate 
on this more generic feature, (rather than on more exotic as- 




Hole concencentration x 

FIG. 5: Typical phase diagram of hole-doped high T c superconduc- 
tors. There exists a pseudogap phase above T c in the underdoped 
regime. Here SC denotes superconductor, and T* is the temperature 
at which the pseudogap smoothly turns on. 
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FIG. 6: Temperature dependence of the excitation gap at the antin- 
odal point (7r,0) in Bi 2 Sr 2 CaCu208+i (BSCCO) for three different 
doping concentrations from near-optimal (discs) to heavy underdop- 
ing (inverted triangles), as measured by angle-resolved photoemis- 
sion spectroscopy (from Ref.|3J). 



pects), which they have in common with other superconduc- 
tors in their distinctive class. 

In Fig. [5] we show a sketch of the phase diagram for the 
hole-doped copper oxide superconductors. Here x represents 
the concentration of holes which can be controlled by, say, 
adding Sr substitutionally to Lai_ a; Sr 2 ;Cu04. At zero and 
small x the system is an antiferromagnetic (AFM) insulator. 
Precisely at half filling (x = 0) we understand this insulator 
to derive from Mott effects. These Mott effects may or may 
not be the source of the other exotic phases indicated in the 
diagram, i.e., the superconducting (SC) and the "pseudogap" 
(PG) phases. Once AFM order disappears the system remains 
insulating until a critical hole concentration (typically around 
a few percent) when an insulator-superconductor transition is 
encountered. Here photoemission studies yl |4J] suggest that 
once this line is crossed, fi appears to be positive. For x < 0.2, 
the superconducting phase has a non-Fermi liquid (or pseudo- 
gapped) normal state |2]. We note an important aspect of this 
phase diagram at low x. As the pseudogap becomes stronger 
(T* increases), superconductivity as reflected in the magni- 
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FIG. 7: Temperature dependence of fermionic excitation gaps A and superfluid density p a for various doping concentrations (from Ref. E^h . 
When A(T ) 7^ 0, there is little correlation between A(T) and p s (T); this figure suggests that something other than fermionic quasi-particles 
(e.g., bosonic excitations) may be responsible for the disappearance of superconductivity with increasing T. Figure (b) shows a quasi-universal 
behavior for the slope dp s /dT at different doping concentrations, despite the highly non-universal behavior for A(T). 



tude of T c becomes weaker. 

Figure [6] indicates the temperature dependence of the ex- 
citation gap for three different hole stoichiome tries. These 
data 01 were taken from angle resolved photoemission spec- 
troscopy (ARPES) measurement. For one sample shown as 
circles, (corresponding roughly to "optimal" doping) the gap 
vanishes roughly at T c as might be expected for a BCS su- 
perconductor. At the other extreme are the data indicated by 
inverted triangles in which an excitation gap appears to be 
present up to room temperature, with very little temperature 
dependence. This is what is referred to as a highly under- 
doped sample (small x), which from the phase diagram can 
be seen to have a rather low T c . Moreover, T c is not evident 
in these data on underdoped samples. 

While the high T c community has focused on pseudogap 
effects above T c , there is a good case to be made that these ef- 
fects also persist below. STM data 1 33] taken below T c within 
a vortex core indicate that there is a clear depletion of the den- 
sity of states around the Fermi energy in the normal phase 
within the core. These data underline the fact that the exis- 
tence of an energy gap has little or nothing to do with the ex- 
istence of phase coherent superconductivity. It also underlines 
the fact that pseudogap effects effectively persist below T c ; the 
normal phase underlying superconductivity for T < T c is not 
a Fermi liquid. 

Analysis of thermodynamical data 14311 has led to a sim- 
ilar inference. For the PG case, the entropy extrapolated into 
the superfluid phase, based on Fermi liquid theory, becomes 
negative. In this way Loram and co-workers |43] deduced 
that the normal phase underlying the superconducting state 
is not a Fermi liquid. Rather, they claimed to obtain proper 
thermodynamics, it must be assumed that this state contains 
a persistent pseudogap. In this way they argued for a dis- 
tinction between the excitation gap A and the superconduct- 
ing order parameter, within the superconducting phase. To 
fit their data they presume a modified fermionic dispersion 
E k = v/(e k -/i) 2 + A 2 (T) where 

A 2 (T) = A 2 C (T) + A 2 pg (4) 

Here A pg is taken on phenomenological grounds to be T- 
independent. While Eq. is also found in BCS-BEC 



crossover theory, there are important differences. In the lat- 
ter approach A pg — > as T — > 0. 

Finally, Fig.Qmakes the claim for a persistent pseudogap 
below T c in an even more suggestive way. Figure [7^ repre- 
sents a schematic plot of excitation gap data such as are shown 
in Fig.|6] Here the focus is on temperatures below T c . Most 
importantly, this figure indicates that the T dependence in A 
varies dramatically as the stoichiometry changes. Thus, in the 
extreme underdoped regime, where PG effects are most in- 
tense, there is very little T dependence in A below T c . By 
contrast at high x, when PG effects are less important, the 
behavior of A follows that of BCS theory. What is most im- 
pressive however, is that these wide variations in A(T) are 
not reflected in the superfluid density p s (T). Figure then 
indicates that, despite the highly non-universal behavior for 
A(T), the superfluid density does not make large excursions 
from its BCS- predicted form. This is difficult to understand if 
the fermionic degrees of freedom through A(T) are dominat- 
ing at all x. Rather this figure suggests that something other 
than fermionic excitations is responsible for the disappearance 
of superconductivity, particularly in the regime where A(T) 
is relatively constant in T. At the very least pseudogap effects 
must persist below T c . 

The phase diagram also suggests that pseudogap effects be- 
come stronger with underdoping. How does one accommo- 
date this in the BCS-BEC crossover scenario? At the simplest 
level one may argue that as the system approaches the Mott 
insulating limit, fermions are less mobile and the effective- 
ness of the attraction increases. In making the connection be- 
tween the strength of the attraction and the variable x in the 
cuprate phase diagram we will argue that it is appropriate to 
simply fit T*(x). In this Review we do not emphasize Mott 
physics because it is not particularly relevant to the atomic 
physics problem. It also seems to be complementary to the 
BCS-BEC crossover scenario. It is understood that both com- 
ponents are important in high T c superconductivity. It should 
be stressed that hole concentration x in the cuprates plays the 
role of applied magnetic field in the cold atom system. These 
are the external parameters which serve to tune the BCS-BEC 
crossover. 

Is there any evidence for bosonic degrees of freedom in 



the normal state of high T c superconductors? The answer 
is unequivocally yes: meta-stable bosons are observable as 
superconducting fluctuations. These effects are enhanced in 
the presence of the quasi-two dimensional lattice structure of 
these materials. In the underdoped case, one can think of T* 
as marking the onset of preformed pairs which are closely 
related to fluctuations of conventional superconductivity the- 
ory, but which are made more robust as a result of BCS-BEC 
crossover effects, that is, stronger pairing attraction. A num- 
ber of people have argued 1451 14611 that fluctuating normal 
state vortices are responsible for the anomalous transport be- 
havior of the pseudogap regime. It has been proposed 14711 that 
these data may alternatively be interpreted as suggesting that 
bosonic degrees of freedom are present in the normal state. 



D. Summary of Cold Atom Experiments: Crossover in the 
Presence of Feshbach Resonances 

There has been an exciting string of developments over the 
past few years in studies of ultracold fermionic atoms, in par- 
ticular, 6 Li and 40 K, which have been trapped and cooled 
via magnetic and optical means. Typically these traps con- 
tain 10 5 atoms at very low densities w 10 13 cm~ 3 . Here 
the Fermi temperature in a trap can be estimated to be of 
the order of a micro-Kelvin. It was argued on the basis of 
BCS theory alone 13, and rather early on (1997), that the 
temperatures associated with the superfluid phases may be 
attainable in these trapped gases. This set off a search for 
the "holy grail" of fermionic superfluidity. That a Fermi de- 
generate state could be reached at all is itself quite remark- 
able; this was was first reported 14 911 by Jin and deMarco in 
1999. By late 2002 reports of unusual hydrodynamics in a 
degenerate Fermi gas indicated that strong interactions were 
present 1 50] . This strongly interacting Fermi gas (associated 
with the unitary scattering regime) has attracted widespread 
attention independent of the search for superfluidity, because 
it appears to be a prototype for analogous systems in nuclear 
physics IdII 15211 and in quark-gluon plasmas 153115411 . More- 
over, there has been a fairly extensive body of analytic work 
on the ground state properties of this regime B55l 15611 . which 
goes beyond the simple mean field wave function ansatz. 

As a consequence of attractive s-wave interactions between 
fermionic atoms in different hyperfine states, it was antici- 
pated that dimers could also be made. Indeed, these molecules 
formed rather efficiently |58, 59, 60] as reported in mid-2003 
either via three body recombination 1 6 1 ] or by sweeping the 
magnetic field across a Feshbach resonance. Moreover, they 
are extremely long lived [59]. From this work it was rel- 
atively straightforward to anticipate that a Bose condensate 
would also be achieved. Credit goes to theorists such as Hol- 
land and to Griffin 1 29] and their co-workers for recog- 
nizing that the superfluidity need not be only associated with 
condensation of long lived bosons, but in fact could also de- 
rive, as in BCS, from fermion pairs. In this way, it was argued 
that a suitable tuning of the attractive interaction via Feshbach 
resonance effects, would lead to a realization of a BCS-BEC 
crossover. 
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FIG. 8: (Color) Spatial density profiles of a molecular cloud of 
trapped 40 K atoms in the BEC regime in the transverse directions 
after 20 ms of free expansion (from Ref. 57), showing thermal molec- 
ular cloud above T c (left) and a molecular condensate (right) below 
T c . (a) shows the surface plots, and (b) shows the cross-sections 
through images (dots) with bimodal fits (lines). 

By late 2003 to early 2004, four groups j53E3|6l|6l had 
observed the "condensation of weakly bound molecules" (that 
is, on the a s > side of resonance), and shortly thereafter a 
number had also reported evidence for superfluidity on the 
BCS side (6^, EE E3, E3] . The BEC side is the more straight- 
forward since the presence of the superfluid is reflected in a 
bi-modal distribution in the density profile. This is shown in 
Fig.|8]from Ref. 57, and is conceptually similar to the behav- 
ior for condensed Bose atoms [69]. On the BEC side but near 
resonance, the estimated T c is about O.STp, with condensate 
fractions varying from 20% or so to nearly 100%. The con- 
densate lifetimes are relatively long in the vicinity of reso- 
nance, and fall off rapidly as one goes deeper into the BEC. 
However, for a s < there is no clear expectation that the den- 
sity profile will provide a signature of the superfluid phase. 

These claims that superfluidity may have been achieved on 
the BCS side {a s < 0) of resonance were viewed as particu- 
larly exciting. The atomic community, for the most part, felt 
the previous counterpart observations on the BEC side were 
expected and not significantly different from condensation in 
Bose atoms. The evidence for this new form of "fermionic su- 
perfluidity" rests on studies I65tl66ll that perform fast sweeps 
from negative a s to positive a s across the resonance. The field 
sweeps allow, in principle, a pairwise projection of fermionic 
atoms (on the BCS side) onto molecules (on the BEC side). It 
is presumed that in this way one measures the momentum dis- 
tribution of fermion pairs. The existence of a condensate was 
thus inferred. Other experiments which sweep across the Fes- 
hbach resonance adiabatically, measure the size of the cloud 
after release lr34ll or within a trap l70ll . 

Evidence for superfluidity on the BCS side, which does not 
rely on the sweep experiments, has also been deduced from 
collective excitations of a fermionic gas 1671 17 ill. Pairing gap 
measurements with radio frequency (RF) spectroscopy probes 
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1 68] have similarly been interpreted as providing support 
for the existence of superfluidity, although more directly these 
experiments establish the existence of fermion pairs. Quite 
recently, evidence for a phase transition has been presented 
via thermodynamic measurements and accompanying theory 
0. The latter, like the theory iQ of RF experiments El, 
is based on the formalism presented in this Review. A most 
exciting and even more recent development has been the ob- 
servation of vortices 1 74] which appears to provide a smoking 
gun for the existence of the superfiuid phase. 



II. THEORETICAL FORMALISM FOR BCS-BEC 
CROSSOVER 

A. Many-body Hamiltonian and Two-body Scattering Theory 

We introduce the Hamiltonian 1II7I I29I l75ll used in the 
cold atom and high T c crossover studies. The most general 
form for this Hamiltonian consists of two types of interac- 
tion effects: those associated with the direct interaction be- 
tween fermions parametrized by U, and those associated with 
"fermion-boson" interactions, whose strength is governed by 



k,<r q 

+ X C/ ( k ' k ')« q /2+k : t (l q/2-k4 a q/2-kM a q/2+k',T 
q,k,k' 

+ J2 (3( k ) 6 q a q/2-k,i a q/2+k,T + h ' C ) ^ 



q,k 



Here the fermion and boson kinetic energies are given by — 
k 2 /2m, and e™ b — q 2 /2M, and v is an important parameter 
which represents the magnetic field-induced "detuning". Here 
we use the convention h = ks = c = 1. In this two channel 
problem the ground state wavefunction is slightly modified 
and given by 



*0 = *0 ® *0 



(6) 



where the molecular or Feshbach boson contribution is as 
given in Ref.l76l 

Whether both forms of interactions are needed in either sys- 
tem is still under debate. The bosons (6^) of the cold atom 
problem iITtL fT8il are referred to as belonging to the "closed 
channel". These spin-singlet molecules represent a separate 
species, not to be confused with the ("open channel") fermion 
pairs (a^a^kX which are associated with spin triplet. As a 
result of virtual occupation of the bound state of the closed 
channel the interaction between open channel fermions can 
be tuned (through applied magnetic field) to vary from weak 
to very strong. 

In this review we will discuss the behavior of crossover 
physics both with and without the closed-channel. Previous 
studies of high T c superconductors have invoked a similar 
bosonic term 127L l77l l78l l79ll as well, although less is known 
about its microscopic origin. This fermion-boson coupling is 



not to be confused with the coupling between fermions and 
a "pairing-mechanism"-related boson {[b + W]a)a) such as 
phonons in a metal superconductor. The coupling b^aa and 
its Hermitian conjugate represent a form of sink and source 
for creating fermion pairs, in this way inducing superconduc- 
tivity in some ways, as a by-product of Bose condensation. 

It is useful at this stage to introduce the s-wave scattering 
length, a, defined by the low energy limit of two body scatter- 
ing in vacuum. We begin with the effects of U only, presum- 
ing that U is always an attractive interaction (U < 0) which 
can be arbitrarily varied, 



m 
4:T:a 



1 x 1 



2e k 



(7) 



We may define a critical value U c of the potential as that asso- 
ciated with the binding of a two particle state in vacuum. We 
can write down an equation for U c given by 



u- 1 = 



E- 



(8) 



although specific evaluation of U c requires that there be a 
cut-off imposed on the above summation, associated with 
the range of the potential. The fundamental postulate of 
crossover theory is that even though the two-body scatter- 
ing length changes abruptly at the unitary scattering condi- 
tion (\a\ = 00), in the N-body problem the superconductivity 
varies smoothly through this point. 

Provided we redefine the appropriate "two body" scattering 
length, Equation Q holds even in the presence of Feshbach 
effects 1 28, 29]. It has been shown that U in the above equa- 
tions is replaced by 



U 



Ueff^U 



9 



2/j, — v 



(9) 



and we write a — ► a*. Experimentally, the two body scattering 
length a* varies with magnetic field B. Thus we have 



1 



Ana* 



U, 



eff 



y- 



2e k 



(10) 



More precisely the effective interaction between two fermions 
is momentum and energy dependent. It arises from a second 
order process involving emission and absorption of a closed- 
channel molecular boson. The net effect of the direct plus 
indirect interactions is given by U e tf(Q) = U + g 2 Do(Q), 
where Dq(Q) = l/[iVL n — e™ h — v + 2n] is the non-interacting 
molecular boson propagator. Here and throughout we use 
a four-momentum notation, Q = (q, zf2„), and its analyt- 
ical continuation, Q — > (q, ft + i0 + ), and write = 
Tj^n Sq' where Sl„ is a Matsubara frequency. What ap- 
pears in the gap equation, however, is U e ff(Q — 0) which 
we define to be U e ff. When the open-channel attraction U 
is weak, clearly, 2p < v is required so that the Feshbach- 
induced interaction is attractive. In the extreme BEC limit 
v = 2/i. However, when a deep bound state exists in the 
open channel, such as in 40 K, the system may evolve into a 
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FIG. 9: Characteristic behavior of the scattering length for 6 Li in the 
three regimes. 



metastable state such that 2/i > v in the BEC regime and 
there is a point on the BCS side where U e f / = precisely. 

Figure [5] presents a plot of this scattering length kpa s = 
kpa* for the case of 6 Li. It follows that a s is negative when 
there is no bound state, it tends to — oo at the onset of the 
bound state and to +oo just as the bound state stabilizes. It 
remains positive but decreases in value as the interaction be- 
comes increasingly strong. The magnitude of a s is arbitrarily 
small in both the extreme BEC and BCS limits, but with op- 
posite sign. 



B. T- Matrix-Based Approaches to BCS-BEC Crossover in the 
Absence of Feshbach Effects 



To address finite temperature in a way which is consistent 
with Eq. Q, or with alternative ground states, one introduces 
a T-matrix approach. Here one solves self consistently for 
the single fermion propagator (G) and the pair propagator (i). 
That one stops at this level without introducing higher order 
Green's functions (involving three, and four particles, etc) is 
believed to be adequate for addressing a leading order mean 
field theory such as that represented by Eq. Q. One can see 
that pair-pair (boson-boson) interactions are only treated in 
a (generalized) mean field averaging procedure; they arise ex- 
clusively from the fermions and are sufficiently weak so as not 
to lead to any incomplete condensation in the ground state, as 
is compatible with Eq. Q. 

In this section we demonstrate that at the T-matrix level 
there are three distinct schemes which can be implemented to 
address BCS-BEC crossover physics. Above T c , quite gener- 
ally one writes for the i-matrix 



t(Q) 



u 



i + u x (Q) 



(ii) 



and theories differ only on what is the nature of the pah- 
susceptibility x(Q)> an( l tne associated self energy of the 
fermions. Below T c one can also consider a T-matrix ap- 
proach to describe the particles and pairs in the condensate. 
For the most part we will defer extensions to the broken sym- 
metry phase to Section liTcl 



In analogy with Gaussian fluctuations, Nozieres and 
Schmitt-Rink considered 1 24] 



xo(Q) = J2 G o( K ) G o(Q-K) 



(12) 



A" 



with self energy 

Z (K) = J2t(Q)Go(Q ~ K) , (13) 



where Go(K) is the noninteracting fermion Green's function. 
The number equation of the Nozieres Schmitt-Rink scheme 
nam is then deduced in an approximate fashion 1 80 ] by 
using a leading order series for G with 



G — Go + GoSqGo 



(14) 



It is straightforward, however, to avoid this approximation in 
Dyson's equation, and a number of groups | JjJ 13711 have ex- 
tended NSR in this way. 
Similarly one can consider 



X 



(Q) = J2G(K)G(Q-K) 



K 



with self energy 



Q 



(15) 



(16) 



This latter scheme has been also extensively discussed in the 
literature, by among others, Haussmann |81], Tchernyshyov 
l82ll and Yamada and Yanatse l4flll . 

Finally, we can contemplate the asymmetric form l25ll for 
the T-matrix, so that the coupled equations for t(Q) and 
G(K) are based on 



x(Q)=Y, G ( K ) G o(Q-K) 



K 



with self energy 



£(*0 = X>(Q)G (Q - *Q 



(17) 



(18) 



It should be noted, however, that this asymmetric form can be 
derived from the equations of motion by truncating the infi- 
nite series at the three particle level, G3, and then factorizing 
the G3 into one- and two-particle Green's functions 18311 . The 
other two schemes are constructed diagrammatically or from a 
generating functional, (as apposed to derived from the Hamil- 
tonian). 

It will be made clear in what follows that, if one's goal is 
to extend the usual crossover ground state of Eq. to finite 
temperatures, then one must choose the asymmetric form for 
the pair susceptibility, as shown in Eq. M7\ . Other approaches 
such as the NSR approach to T c , or that of Haussmann lead 
to different ground states which should, however, be very in- 
teresting in their own right. These will need to be charac- 
terized in future. Indeed, the work of Strinati group has also 
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emphasized that the ground state associated with the T c calcu- 
lations based on NSR is distinct from that in the simple mean 
field theory of Eq. (HJ, and they presented some aspects of this 
comparison in Ref.l84l 

Other support for this GGo -based T-matrix scheme comes 
from its equivalence to self consistent Hartree-approximated 
Ginzburg-Landau theory 1 85 ] . Moreover, there have been de- 
tailed studies 1 86] to demonstrate how the superfluid density 
p s can be computed in a fully gauge invariant (Ward Identity 
consistent) fashion, so that it vanishes at the self consistently 
determined T c . Such studies are currently missing for the case 
of the other two T-matrix schemes. 



C. Extending conventional Crossover Ground State to T ^ 0: 
T-matrix scheme in the presence of closed-channel molecules 

In the T-matrix scheme we employ, the pairs are de- 
scribed by the pair susceptibility x(Q) = J2k^o(Q — 
K)G(K)ip%._, 2 where G depends on a BCS-like self en- 
ergy S(JsT) w -A 2 Go(-K)ip1. Throughout this section 
<£>k = cxp{ — k 2 /2kfy introduces a momentum cutoff, where 
fco represents the inverse range of interaction, which is as- 
sumed infinite for a contact interaction. 

The noncondensed pairs 18711 have propagator t pg (Q) = 
U e ff(Q)/[l + U e ff(Q)x{Q)}, where U e ff is the effective 
pairing interaction which involves the direct two-body interac- 
tion U as well as virtual excitation processes associated with 
the Feshbach resonance ll29l. l87ll . At small Q, t pg can be ex- 
panded, after analytical continuation (iQ n — > 17 + i0 + ), as 



the first of several constraints needed to close the set of equa- 
tions, (i) One has a constraint on the total number of pairs Il75ll 
which can be viewed as analogous to the usual BEC number 
constraint 



A 2 (T) = A 2 C (T) + A 2 AT) . 



(20) 



To determine A, (ii) one imposes the BEC-like constraint 
that the pair chemical potential vanishes in the superfluid 
phase: 



T <T C 



This yields 



so that 



t p - g 1 (Q^Q) = Q = U- f 1 f (0) + x(0) 



U effW + 1^ 2K ^ k = ° 



(21) 



(22) 



(23) 



Importantly, below T c , A satisfies the usual BCS gap equa- 
ti on. Here we introd uce the quasiparticle dispersion = 

V (ek — m) 2 + A 2( Pk> anc ^ f( x ) * s tne F errm distribution 
function. 

(iii) In analogy with the standard derivation of BEC, the 
total contribution of noncondensed pairs is readily computed 
by simply adding up their number, based on the associated 
propagator 



a p 9 = -Ewq) 



(24) 



tpg{Q) 



n-n 



q T" Impair 



(19) 



The parameters appearing in Eq. \\9\ are discussed in more 
detail in Ref. HI Here Z' 1 is a residue and Q q = q 2 /2M* 
the pair dispersion, where M* is the effective pair mass. The 
latter parameter as well as the pair chemical potential p va i r 
depends on the important, but unknown, gap parameter A 
through the fermion self energy £, The decay width Tq is 
negligibly small for small Q below T c . 

While there are alternative ways of deriving the self con- 
sistent equations which we use, (such as a decoupling of the 
Green's function equations of motion |83]), here we present 
an approach which shows how this GGo-based T-matrix 
scheme has strong analogies with the standard theory of BEC. 
But, importantly this BEC is embedded in a self consistent 
treatment of the fermions. Physically, one should focus on 
A as reflecting the presence of bosonic degrees of freedom. 
In the fermionic regime (p > 0), it represents the energy re- 
quired to break the pairs, so that A is clearly associated with 
the presence of "bosons". In the bosonic regime, A 2 directly 
measures the density of pairs. 

In analogy with the standard theory of BEC, it is expected 
ll87ll that A contains contributions from both noncondensed 
and condensed pairs. The associated densities are propor- 
tional to A 2 AT) and A 2 C (T), respectively. We may write 



One can rewrite Eq. J24i so that it looks more directly like 
a number equation, by introducing the Bose distribution func- 
tion b(x) for noncondensed pairs as 



K a = z 



Kn g ,T), 



(25) 



so that the noncondensed pair density is given by ZA 2 g . Note 
that the right hand sides of the previous two equations depend 
on the unknown A through the self energy appearing in G 
which, in turn enters t pg or il q . Also note that at T = 0, 
A pg = so that all pairs are condensed as is consistent with 
the mean-field BCS-Leggett ground state. 

Finally, in analogy with the standard derivation of BEC, 
(iv) one can then compute the number of condensed pairs as- 
sociated with A sc , given that one knows the total A and the 
noncondensed component. 

Despite this analogy with BEC, fermions are the fundamen- 
tal particles in the system. It is their chemical potential p that 
is determined from the number conservation constraint 



2n 



b0 



2nh = n 



f 



2nl 



(26) 



Here n;,o an d n b represent the density of condensed and non- 
condensed closed-channel molecules, respectively, njj ot is the 
sum, and nj = 2^G(X) is the atomic density associ- 
ated with the open-channel fermions. These closed channel 
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FIG. 10: Typical behavior of T c as a function of l/ZcfO in a ho- 
mogeneous system. T c follows the BCS predictions and approaches 
the BEC asymptote 0.2182> in the BEC limit. In the intermediate 
regime, it reaches a maximum around l/kpa = and a minimum 
around where [i = 0. 



FIG. 1 1 : Typical behavior of T c of a Fermi gas in a trap as a function 
of l/kpa. It follows BCS prediction in the weak coupling limit, 
l/kpa << —1, and approaches the BEC asymptote 0.518Tf in the 
limit 1/fc.pct — » +oo. In contrast to the homogeneous case in Fie. llQI 
the BEC asymptote is much higher due to a compressed profile for 
trapped bosons. 



fermions have a propagator D(Q), which we do not discuss 
in much detail in order to make the presentation simpler. Here 

n 6 = -^D(Q)«^^6(O q ), (27) 
Q q 

where b(x) is the Bose distribution function. The renor- 
malized propagator D(Q) is given by the same equation as 
Eq. ( I19> with a different residue Z^ 1 — > Z\>. 

In this way the system of equations is complete 1 87 ] . The 
numerical scheme is straightforward in principle. We compute 
A (and /i) via Eqs. d23i and (I26> . to determine the contribution 
from the condensate A sc via Eqs. d20l and (1241 . Above T c , 
this theory must be generalized to solve self consistently for 
Impair which no longer vanishes ll88ll . 

III. PHYSICAL IMPLICATIONS: ULTRACOLD ATOM 
SUPERFLUIDITY 

In this section we compare four distinct classes of exper- 
iments on ultracold trapped fermions with theory. These are 
thermodynamics [73, 88], temperature dependent density pro- 
files J89ll . RF pairing gap spectroscopy [68, 90l l9lll . and col- 
lective mode measurements 11671 17 ill . We address all four 
experiments in the context of the mean field ground state of 
Eq. Q, and its finite temperature extension discussed in Sec- 
tion III CI That there appears to be good agreement between 
theory and experiment lends rather strong support to the sim- 
ple mean field theory, which is at the center of this Review. 
Interestingly, pseudogap effects are evident in various ways in 
these experiments and this serves to tie the ultracold fermions 
to the high T c superconductors. 

A. T c Calculations and Trap Effects 

Before turning to experiment, it is important to discuss the 
behavior of the transition temperature which is plotted as a 



function of scattering length in Fig.[l0|for the homogeneous 
case, presuming s-wave pairing. We discuss the effects of d- 
wave pairing in Section ffvl in the context of application to 
the cuprates. Starting from the BCS regime this figure shows 
that T c initially increases as the interaction strength increases. 
However, this increase competes with the opening of a pseu- 
dogap or excitation gap A(T C ). Technically, the pairs become 
effectively heavier before they form true bound states. Even- 
tually T c reaches a maximum (very near unitarity) and then 
decreases slightly until field strengths corresponding to the 
point where /i becomes zero. At this field value (essentially 
where T c is minimum), the system becomes a "bosonic" su- 
perfluid, and beyond this point T c increases slightly to reach 
the asymptote corresponding to an ideal Bose gas. Trap 
effects change these results only quantitatively as seen in 
Fig. ^2 To treat these trap effects one introduces the local 
density approximation (LDA) in which T c is computed under 
the presumption that the chemical potential /.i — ► fj, — V(r) 
. Here we consider a spherical trap with V(r) = ^mui 2 r 2 . 
The Fermi energy Ep is determined by the total atom num- 
ber N via E F = k B T F = huj^N) 1 ^ 3 = h 2 k 2 F /2m, where 
hp is the Fermi wavevector at the center of the trap. It can 
be seen that the homogeneous curve is effectively multiplied 
by an "envelope" curve when a trap is present. This envelope, 
with a higher BEC asymptote, reflects the fact that the particle 
density at the center of the trap is higher in the bosonic, rela- 
tive to the fermionic case. In this way T c is relatively higher 
in the BEC regime, as compared to BCS, whenever a trap is 
present. 

Figure EJpresents a plot of the position dependent excita- 
tion gap A(r) and particle density n{r) profile over the extent 
of the trap. An important point needs to be made: because 
the gap is largest at the center of the trap, bosonic excitations 
will be dominant there. At the edge of the trap, by contrast, 
where fermions are only weakly bound (since A(r) is small), 
the excitations will be primarily fermionic. We will see the 
implications of these observations as we examine thermody- 
namic and RF spectra data in the ultracold gases. 
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FIG. 12: Typical spatial profile of T — density n(r) and fermionic 
excitation gap A(r) of a Fermi gas in a trap. The curves are com- 
puted at unitarity, where 1/fcpci = 0. Here Rtf is the Thomas- 
Fermi radius. 



B. Thermodynamical Experiments 

Figure [O] present a plot which compares experiment and 
theory in the context of thermodynamic experiments f73l l88il 
on trapped fermions. Plotted on the vertical axis is the en- 
ergy which can be input in a controlled fashion experimen- 
tally. The horizontal axis is temperature which is calibrated 
theoretically based on an effective temperature T introduced 
phenomenologically, and discussed below. The experimen- 
tal data are shown for the (effectively) non-interacting case as 
well as unitary. In this discussion we treat the non-interacting 
and BCS cases as essentially equivalent since A is so small 
on the scale of the temperatures considered. The solid curves 
correspond to theory for the two cases. Although not shown 
here, even without a temperature calibration, the data suggests 
a phase transition is present in the unitary case. This can be 
seen as a result of the change in slope of E(T) as a function 
of f. 

The phenomenological temperature T is relatively easy to 
understand. What was done experimentally to deduce this 
temperature was to treat the unitary case as an essentially free 
Fermi gas to, thereby, infer the temperature from the width 
of the density profiles, but with one important proviso: a nu- 
merical constant is introduced to account for the fact that the 
density profiles become progressively narrower as the system 
varies from BCS to BEC. This systematic variation in the pro- 
file widths reflects the fact that in the free Fermi gas case, 
Pauli principle repulsion leads to a larger spread in the particle 
density than in the bosonic case. And the unitary regime has 
a profile width which is somewhere in between, so that one 
parametrizes this width by a simple function of ft. We can 
think of ft as reflecting bosonic degrees of freedom, within 
an otherwise fermionic system. At ft = the system is a 
free Fermi gas. The principle underlying this rescaling of the 
non-interacting gas is known as the "universality hypothesis" 
ll92l.l93ll . At unitarity, the Fermi energy of the non-interacting 
system is the only energy scale in the problem (for the widely 
used contact potential) since all other scales associated with 
the two-body potential drop out when a s — > ±00. We re- 
fer to this phenomenological fitting temperature procedure as 
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FIG. 13: (color) Energy E vs physical temperature T. The upper 
curve and data points correspond to the BCS or essentially free Fermi 
gas case, and the lower curve and data correspond to unitarity. The 
latter provide indications for a phase transition. The inset shows how 
temperature must be recalibrated below T c . From Ref. l73ll . 



Thomas-Fermi (TF) fits. 

An interesting challenge was to relate this phenomenolog- 
ical temperature T to the physical temperature T; more pre- 
cisely one compares \/l + ftT and T. This relationship is 
demonstrated in the inset of Fig.^] And it was in this way 
that the theory and experiment could be plotted on the same 
figure, as shown in the main body of Fig.^] The inset was 
obtained using theory only. The theoretically produced pro- 
files were analyzed just as the experimental ones to extract 
\J\ + ftT and compare it to the actual T, Above T c no re- 
calibration was needed as shown by the straight line going 
through the diagonal. Below T c the phenomenologically de- 
duced temperatures were consistently lower than the physical 
temperature. That the normal state temperatures needed no 
adjustment shows that the phenomenology captures important 
physics. It misses, however, an effect associated with the pres- 
ence of a condensate which we will discuss shortly. 

We next turn to a more detailed comparison of theory and 
experiment for the global and low T thermodynamics. Fig- 
ure^]presents a blow-up of E at the lowest T comparing the 
unitary and non-interacting regimes. The agreement between 
theory and experiment is quite good. In the figure, the tem- 
perature dependence of E reflects primarily fermionic excita- 
tions at the edge of the trap, although there is a small bosonic 
contribution as well. It should be noted that the theoretical 
plots were based on fitting ft to experiment by picking a mag- 
netic field very slightly off resonance. [In the simple mean 
field theory ft = —0.41, and in Monte Carlo simulations i5^1 
ft = —0.54. Both these theoretical numbers lie on either side 
of experiment 1 73] where ft — —0.49]. 

Figure ^] presents a wider temperature scale plot which, 
again, shows very good agreement. Importantly one can see 
the effect of a pseudogap in the unitary case. The tempera- 
ture T* can be picked out from the plots as that at which the 
non-interacting and unitary curves intersect. This corresponds 
roughly to T* w 2T C . 
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FIG. 14: (color) Low temperature comparison of theory (curves) and 
experiments (symbols) in terms of E/Ef (Ef = fesTp) per atom 
as a function of T/Tf, for both unitary and noninteracting gases in 
a Gaussian trap. From Ref.Q. 
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FIG. 15: (color) Same as Fig. ll4l but for a much larger range of tem- 
perature. The quantitative agreement between theory and experiment 
is very good.The fact that the two experimental (and the two theoret- 
ical) curves do not merge until higher T* > T c is consistent with the 
presence of a pseudogap. 



C. Temperature Dependent Particle Density Profiles 

In order to understand more deeply the behavior of the ther- 
modynamics, we turn next to a comparison of finite T density 
profiles. Experiments which measure these profiles [7(1 19411 
all report that they are quite smooth at unitarity, without any 
signs of the bimodality seen in the BEC regime. We discuss 
these profiles in terms of the four panels in Fig.^] These fig- 
ures are a first step in understanding the previous temperature 
calibration procedure. 

In this figure we compare theory and experiment for the 
unitary case. The experimental data were estimated to corre- 
spond to roughly this same temperature (T/Tp = 0.19) based 
on the calibration procedure discussed above. The profiles 
shown are well within the superfluid phase (T c w 0.3Tp at 
unitarity). This figure presents Thomas-Fermi fits [94] to (a) 
the experimental and (b) theoretical profiles as well as (c) their 
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FIG. 16: (color) Temperature dependence of (a) experimental one- 
dimensional spatial profiles (circles) and TF fit (line) from Ref . l94l 

(b) TF fits (line) to theory both at T « 0.7T C » 0.19T F (circles) and 

(c) overlay of experimental (circles) and theoretical (line) profiles, as 
well as (d) relative rms deviations (\ 2 ) associated with these fits to 
theory at unitarity. The circles in (b) are shown as the line in (c). The 
profiles have been normalized so that N = f n(x)dx = 1, and we 
set Rtf = 100 /J,m in order to overlay the two curves, \ 2 reaches a 
maximum around T = 0.19Tf. 



comparison, for a chosen Rtf = 100 /im, which makes it 
possible to overlay the experimental data (circles) and theo- 
retical curve (line). Finally Fig.ll6H indicates the relative x 2 
or root-mean-square (rms) deviations for these TF fits to the- 
ory. This figure was made in collaboration with the authors 
of Ref. 1 94]. Two of the three dimensions of the theoretical 
trap profiles were integrated out to obtain a one-dimensional 
representation of the density distribution along the transverse 
direction: n(x) = J dydzn(r). 

This figure is in contrast to earlier theoretical studies which 
predict a significant kink at the condensate edg e which ap- 
pears not to have been seen experimentally 17(1 19411 . More- 
over, the curves behave monotonically with both temperature 
and radius. Indeed, in the unitary regime the generalized TF 
fitting procedure of Ref. 94 works surprisingly well. And 
these reasonable TF fits apply to essentially all temperatures 
investigated experimentally 19411 . as well as theoretically, in- 
cluding in the normal state. 

It is important to establish why the profiles are so smooth, 
and the condensate is, in some sense, rather invisible, except 
for its effect on the TF-inferred-temperature. This apparent 
smoothness can be traced to the presence of noncondensed 
pairs of fermions which need to be included in any consistent 
treatment. Indeed, these pairs below T c are a natural counter- 
part of the pairs above T c which give rise to pseudogap effects. 

To see how the various contributions enter into the trap pro- 
file, in Fig.^l we pl°t a decomposition of this profile for var- 
ious temperatures from below to above T c . The various color 
codes indicate the condensate along with the noncondensed 
pairs and the fermions. This decomposition is based on the 
superfluid density so that all atoms participate in the conden- 
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FIG. 17: (color) Decomposition of density profiles at various temper- 
atures at unitarity. Here green (light gray) refers to the condensate, 
red (dark gray) to the noncondensed pairs and blue (black) to the ex- 
cited fermionic states. T c — 0.27Tf, and Rtf is the Thomas-Fermi 
radius. 



sation at T = 0. This, then, forms the basis for addressing 
both thermodynamics and RF pairing-gap spectroscopy in this 
Review. 

The figure shows that by T = T c /2 there is a reasonable 
number of excited fermions and bosons. As anticipated earlier 
in Section HTCl the latter are at the trap edge and the former 
in the center. By T = T c the condensate has disappeared and 
the excitations are a mix of fermions (at the edge) and bosons 
towards the center. Indeed, the noncondensed bosons are still 
present by T = 1.5T C , as a manifestation of a pseudogap ef- 
fect. Only for somewhat higher T 2T C do they disappear 
altogether, so that the system becomes a non-interacting Fermi 
gas. 

Two important points should be made. The noncondensed 
pairs clearly are responsible for smoothing out what otherwise 
would be a discontinuity 1 92, 95] between the fermionic and 
condensate contributions. Moreover, the condensate shrinks 
to the center of the trap as T is progressively raised. It is this 
thermal effect which is responsible for the fact that the TF 
fitting procedure for extracting temperature leads to an under- 
estimate as shown in the inset to Fig.^] The presence of the 
condensate tends to make the atomic cloud smaller so that the 
temperature appears to be lower in the TF fits. 



D. RF Pairing gap Spectroscopy 



Measurements |68] of the excitation gap A have been made 
by using a third atomic level, called |3), which does not partic- 
ipate in the superfluid pairing. Under application of RF fields, 
one component of the Cooper pairs, called |2), is presumably 
excited to state |3) . If there is no gap A then the energy it takes 
to excite |2) to |3) is the atomic level splitting uj-23- In the pres- 
ence of pairing (either above or below T c ) an extra energy A 
must be input to excite the state |2), as a result of the breaking 
of the pairs. Figure[^]shows a plot of the spectra near unitar- 
ity for four different temperatures, which we discuss in more 
detail below. In general for this case, as well as for the BCS 
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FIG. 18: Experimental RF Spectra at unitarity. The temperatures la- 
beled in the figure were computed theoretically at unitarity based on 
adiabatic sweeps from BEC. The two top curves, thus, correspond 
to the normal phase, thereby, indicating pseudogap effects. Here 
E F = 2.5^iK, or 52 kHz. From Ref.^l 



and BEC limits, there are two peak structures which appear in 
the data: the sharp peak at CJ23 = which is associated with 
"free" fermions at the trap edge and the broader peak which 
reflects the presence of paired atoms; more directly this broad 
peak derives from the distribution of A in the trap. At high 
T (compared to A), only the sharp feature is present, whereas 
at low T only the broad feature remains. The sharpness of 
the free atom peak can be understood as coming from a large 
phase space contribution associated with the 2^3 excita- 
tions 1 91]. Clearly, these data alone do not directly indicate 
the presence of superfluidity, but rather they provide strong 
evidence for pairing. 

As pointed out in Ref. 1 90] these experiments serve as a 
counterpart to superconducting tunneling in providing infor- 
mation about the excitation gap. A theoretical understanding 
of these data was first presented in Ref. 72 using the frame- 
work of Section ITTCl Subsequent work [91] addressed these 
data in a more quantitative fashion as plotted in Fig.[H)] Here 
the upper and lower panels correspond respectively to inter- 
mediate and low temperatures. For the latter one sees that 
the sharp "free atom" peak has disappeared, so that fermions 
at the edge of the trap are effectively bound at these low 
T. Agreement between theory and experiment is quite sat- 
isfactory, although the total number of particles was adjusted 
somewhat relative to the experimental estimates. 

It is interesting to return to the previous figure (Fig. I18> 
and to discuss the temperatures in the various panels. What 
is measured experimentally are temperatures T' which cor- 
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FIG. 19: Comparison of calculated RF spectra (solid curve, T c ~ 
0.29Tf) with experiment (symbols) in a harmonic trap calculated at 
822 G for the two lower temperatures. The temperatures were chosen 
based on Ref . [ 68] . The particle number was reduced by a factor of 2, 
as found to be necessary in addressing another class of experiments 
(84j. The dashed lines are a guide to the eye. From Ref.l9ll 



respond to the temperature at the start of a sweep from the 
BEC limit to unitarity. Here fits to the BEC-like profiles are 
used to deduce T" from the shape of the Gaussian tails in the 
trap. Based on knowledge about thermodynamics (entropy S), 
and given T', one can then compute the final temperature in 
the unitary regime, assuming S is constant. Indeed, this adi- 
abaticity has been confirmed experimentally in related work 
I70B . We find that the four temperatures are as indicated in the 
figures. Importantly, one can conclude that the first two cases 
correspond to a normal state, albeit close to T c . Importantly, 
these figures suggest that a pseudogap is present as reflected 
by the broad shoulder above the narrow free atom peak. 



E. Collective Breathing Modes at T 

We turn, finally, to a com pari son between theory Ir9rjl 
and experiment 1671 17 U l93l 19811 for the collective breathing 
modes within a trap at T « 0. The very good agreement 
has provided some of the earliest and strongest support for 
the simple mean field theory of Eq. Q. Interestingly, Monte 
Carlo simulations which initially were viewed as a superior 
approach, lead to significant disagreement between theory and 
experiment 1 99] . Shown in Fig.|20]is this comparison for the 
axial mode in the inset and the radial mode in the main body 
of the figure as a function of magnetic field. The experimental 
data are from Ref. 67, The original data on the radial modes 
from Ref. |7l|, was in disagreement with that of Ref. |67j, but 
this has since been corrected 119911 . and there is now a consis- 
tent experimental picture from both the Duke and the Inns- 
bruck groups for the radial mode frequencies. 

At T = 0, calculations of the mode frequencies can be re- 
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FIG. 20: Breathing mode frequencies as a function of k « 
1.695(fc.Fa) -1 , from Tosi et al. 19611 . The main figure and inset plot 
the transverse and axial frequencies, respectively. The solid curves 
are calculations l96ll based on BCS-BEC crossover theory at T = 0, 
and the symbols plot the experimental data from Kinast et al. loTll . 



duced to a calculation of an equation of state for jj, as a func- 
tion of n. One of the most important conclusions from this 
figure is that the behavior in the near-BEC limit (which is still 
far from the BEC asymptote) shows that the mode frequencies 
decrease with incr easin g magnetic field. This is opposite to 
earlier predictions 1 100] based on the behavior of true bosons 
where a Lee-Yang term would lead to an increase. Indeed, the 
pair operators do not obey the commutation relations o f tru e 
bosons except in the zero density or fc^a — > + limit ll 1 ltl . 
Figure [20] thus, underlines the fact that fermionic degrees of 
freedom (or compositeness) are still playing a role at these 
magnetic fields. There are predictions in the literature 1l02ll 
that one needs to achieve fcpa somewhat less than 0.3 (ex- 
perimentally, the smallest values for these experiments are 0.3 
and 0.7 for the various groups) in order to approach the true 
bosonic limit. At this point, then, the simple mean field theory 
will n o longer be adequate. Indeed, there are other indications 
lll03ll of the breakdown of this mean field in the extreme BEC 
limit which are, physically, reflected in the width of the parti- 
cle density profiles. This originates from an overestimate (by 
roughly a factor of 3) of the size of the effective "inter-boson" 
scattering length. 

Overall the mean field theory presented here looks very 
promising. Indeed, the agreement between theory and exper- 
iment is better than one might have anticipated. For the col- 
lective mode frequencies, it appears to be better than Monte 
Carlo calculations lr99Tl . Nevertheless, uncertainties remain. 
Theories which posit a different ground state will need to be 
compared with the four experiments discussed here. It is, fi- 
nally, quite possible that incomplete T = condensation will 
become evident in future experiments. If so, an alte rnative 
wavefunction will have to be contemplated 1 102, 104]. What 
appears to be clear from the current experiments is that, just 
as in high T c superconductors, the ultracold fermionic super- 
fluids exhibit pseudogap effects. These are seen in thermo- 
dynamics, in RF spectra and in the temperature dependence 
of the profiles (through the noncondensed pair contributions). 
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Moreover, while not discussed here, at finite T, damping of 
the collective mode frequencies seems to change qualitatively 
J93ll at a temperature which is close to the estimated T* . 

Looking to the future, at an experimental level, new pairing 
gap spectr oscopies appear to be emerging at a fairly rapid pace 
Il05lll0dl . These will further test the present and subsequent 
theo ries. I ndeed, recently, a probe of the closed channel frac- 
tion 1 106] has been analyzed 1 107] within the present frame- 
work and has led to good quantitative agreement between the- 
ory and experiment. 



IV. PHYSICAL IMPLICATIONS: HIGH T c 
SUPERCONDUCTIVITY 
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A. Phase Diagram and Superconducting Coherence 

The high T c superconductors are different from the ultra- 
cold fermionic superfiuids in one key respect; they are d-wave 
superconductors and their electronic dispersion is associated 
with a quasi-two dimensional tight binding lattice. In many 
ways this is not a profound difference from the perspective of 
BCS-BEC crossover. Figure|^shows a plot of the two impor- 
tant temperatures T c and T* as a function of increasing attrac- 
tive coupling. On the left is BCS and the right is PG. The BEC 
regime is not visible. This is because T c disappears before it 
can be accessed. This disappearance of T c is relatively easy 
to understand. Because the d-wave pairs are more extended 
(than their s-wave counterparts) they experience Pauli princi- 
ple repulsion more intensely. Consequently the pairs localize 
(their mass is infinite) well before the fermionic chemical po- 
tential is negative 1 38 ] . 

The competition between T* and T c , in which as T* in- 
creases, T c decreases, is also apparent in Fig. This is 
a consequence of pseudogap effects. More specifically, the 
pairs become heavier as the gap increases in the fermionic 
spectrum, competing with the increase of T c due to the in- 
creasing pairing strength. It is interesting to compare Fig. 12 1 1 
with the experimental phase diagram plotted as a function of 
the doping concentration x in Fig. [5] If one inverts the hori- 
zontal axis (and ignores the unimportant AFM region) the two 
are very similar. To make an association from coupling U to 
the variable x, it is reasonable to fit T* . It is not particularly 
useful to implement this last step here, since we wish to em- 
phasize crossover effects which are not complicated by "Mott 
physics". 

Because of quasi-two dimensionality, the energy scales of 
the vertical axis in Fig.|^are considerably smaller than their 
three dimensional analogues. Thus, pseudogap effects are 
intensified, just as conventional fluctuation effects are more 
apparent in low dimensional systems. This may be one of 
the reasons why the cuprates are among the first materials to 
clearly reveal pseudogap physics. Moreover, the present cal- 
culations show that in a strictly 2D material, T c is driven to 
zero, by bosonic or fluctuation effects. This is a direct reflec- 
tion of the fact that there is no Bose condensation in 2D. 

The presence of pseudogap effects raises an interesting set 
of issues surrounding the signatures of the transition which 



FIG. 21: Typical phase diagram for a quasi-two dimensional d-wave 
superconductor on a tight-binding lattice at high filling n ~ 0.85 per 
unit cell; here the horizontal axis corresponds to —U/At, where t is 
the in-plane hopping matrix element. 

the high T c community has wrestled with, much as the cold 
atom community is doing today. For a charged superconduc- 
tor there is no difficulty in measuring the superfluid density, 
through the electrodynamic response. Thus one knows with 
certainty where T c is. Nevertheless, people have been con- 
cerned about precisely how the onset of phase coherence is 
reflected in thermodynamics, such as CV or in the fermionic 
spectral function, given that a gap is already present at the 
onset of superconductivity. One understands how phase co- 
herence shows up in BCS theory, since the ordered state is 
always accompanied by the appearance of an excitation gap. 

To addre ss th ese coherence effects one has to introduce a 
distinction 1 1 10] between the self energy associated with non- 
condensed and condensed pairs. This distinction is blurred 
by the approximations made in Section lfl"Cl Within this im- 
proved scheme 1 110] superconducting coherence effects can 
be probed as, presented in Fig. |22l along with a compari- 
son to experiment. Shown are the results of specific heat 
and tunneling calculations and their experimental counter- 
parts I2L I33II . The latter measures, effectively, the density of 
fermionic states. Here the label "PG" corresponds to an ex- 
trapolated normal state in which we set the superconducting 
order parameter A sc to zero, but maintain the the total excita- 
tion gap A to be the same as in a phase coherent, supercon- 
ducting state. Agreement between theory and experiment is 
satisfactory. 

B. Electrodynamics in the superconducting phase 

In some ways the subtleties of phase coherent pairing are 
even more perplexing in the context of electrodynamics. Fig- 
ure Q presents a paradox in which the excitation gap for 
fermions appears to have little to do with the behavior of 
the superfluid density. This superfluid density can be readily 
computed within the BCS-BEC crossover scenario 1251 14411 . 
Particularly important is to include all excitations of the con- 
densate in a fully consistent fashion, compatible with ther- 
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FIG. 22: Extrapolated normal state (PG) and superconducting state (SC) contributions to SIN tunneling and thermodynamics (left), as well as 
comparison with experiments (right) on tunneling for BSCCO 1 33] and on specific heat for Yo.8Cao.2Ba2Cu307_5 1108 1. The theoretical SIN 
curve is calculated for T = T c /2, while the experimental curves are measured outside (dashed line) and inside (solid line) a vortex core. 



modynamics, and which is also manifestly gauge invariant. 
To make contact with electrodynamic experiments, one has to 
introduce the variable x and this is done via a fit to T*(x) 
in the phase diagram. In addition it is also necessary to fit 
p„ (T = 0, x) to experiment, and we do so here, noting that 
1321 the Uemura relation p s ( 0,x) oc T r (x) no longer holds 
for very underdoped samples 1 109, 1 1 ill . By fitting these in- 
dependent quantities we are, in effect accounting for at least 
some aspects of Mott physics. The paradox raised by Fig.0 
is resolved by noting that there are bosonic excitations of the 
condensate |25] and that these become more marked with un- 
derdoping, as pseudogap effects increase. In this way p s does 
not exclusively reflect the fermionic gap, but rather vanishes 
"prematurely" before this gap is zero, as a result of pair exci- 
tations of the condensate. 

This theory can be quantitatively compared with experi- 
ment. Figure [23] presents theoretical and experimental plots 
of the lower critical field, H c i (T), for a group of sev erely un- 
derdoped YBCO crystals as considered in Ref. ll09l There it 
was argued that H C \{T) oc p s (T), so that the lower critical 
field effectively measures the in-plane superfluid density. Ex- 
perimentally what is directly measured is the magnetization 
with applied field parallel to the c-axis. The experimental re- 
sults are shown on the lower two panels and theory on the 
upper two. The left hand figures plot H C %(T) vs T and the 
right hand figures correspond to a rescaling of this function in 
the form H cX (T) /H cl (0) vs T/T c . Theoretically, it is found 
that the fermionic contribution leads to a linear T dependence 
at low T, associated with d-wave pairing, whereas the bosonic 
term introduces a T 3 ' 2 term. Quite remarkably even when the 
Uemura relation no longer holds, there is still a "universality" 
in the normalized plots as shown in both theory and exper- 
iment by the right hand figures. It should be noted that the 



experimental plot contains (at T c = 55.5K) a slightly differ- 
ent cuprate phase known as the ortho-II phase, which does not 
lie on the universal curves. The universality found here can be 
understood as associated with the fact that T c (x), rather than 
A(x), is the fundamental energy scale in p s (T, x). The reason 
that A(x) is not the sole energy scale is that bosonic degrees 
of freedom are also present, and help to drive p K to zer o at 
T c . By contrast, Fermi-liquid based approaches Il09lll 12ll as- 
sume that the fermions are the only relevant excitations, and 
they account for this data by introducing a phenomenological 
parameter a which corresponds to the effective charge of the 
fermionic quasi-particles. 

As anticipated in earlier theoretical calculations l25il the 
bosonic contribution begins to dominate in severely under- 
doped systems so that the slope dH c i / dT (associated with the 
lowest temperatures reached experimentally) should decrease 
with underdoping. Although observed a number of years after 
this prediction, this is precisely what is seen experimentally, 
as shown in Fig. [24] Here the inset plots the experimental 
counterpart data. It can be seen that theory and experiment 
are in reasonably good quantitative agreement. This theoret- 
ical viewpoint is very different from a "Fermi-liquid" based 
treatment of the superconducting state, for which the strong 
decrease in the slope of H c \ was not expected. Within the 
present formalism, the optical conductivity (jl a{u>) is simi- 
larly modified 111 3ll to include bosonic as well as fermionic 
contributions. 



C. Bosonic Power Laws and Pairbreaking Effects 

The existence of noncondensed pair states below T c af- 
fects thermodynamics, in the same way that electrodynamics 
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FIG. 23: (color) Comparison between calculated lower critical field, H c i, as a function of T (upper left panel), and experimental data (lower 
left) from Ref. llfM with variable doping concentration x. The right column shows normalized plots, H c \ (T) / H c \ (0) versus T/T c , for theory 
and experiment, respectively, revealing a quasi-universal behavior with respect to doping, with the exception of the T c = 55. 5K ortho-II 
phase. Both theory plots share the same legends. The quantitative agreement between theory and experiment is quite good. 



is affected, as discussed above. Moreover, one can predict 
ll36ll that the q 2 dispersion will lead to ideal Bose gas power 
laws in thermodynamical and transport properties. These will 
be present in addition to the usual power laws or (for s- 
wave) exponential temperature dependencies associated with 
the fermionic quasi-particles. Note that the q 2 dependence 
is dictated by the ground state of Eq. 0. Clearly this mean 
field like state is inapplicable in the extreme BEC limit, where, 
presumably inter-boson effects become important and lead 
to a linear dispersion. Presumably, in the PG or near-BEC 
regimes, fermionic degrees of freedom are still dominant and 
it is reasonable to apply Eq. Q. Importantly, at present nei- 
ther the cuprates nor the cold atom systems access this true 
BEC regime. 

The consequences of these observations can now be listed 
ll36ll . For a quasi-two dimensional system, C v /T will ap- 
pear roughly constant at the lowest temperatures, although it 
vanishes strictly at T = as T 1 / 2 . The superfluid density 
p s (T) will acquire a T 3 / 2 contribution in addition to the usual 
fermionic terms. By contrast, for spin singlet states, there is 



no explicit pair contribution to the Knight shift. In this way 
the low T Knight shift reflects only the fermions and exhibits 
a scaling with T/A(0) at low temperatures. Experimentally, 
in the cuprates, one usually sees a finite low T contribution to 
C v /T. A Knight shift scaling is seen. Finally, also observed 
is a deviation from the predicted d-wave linear in T power 
law in p s . The new power laws in C v and p s are convention- 
ally attributed to impurity effects. Experiments are not yet 
at a stage to clearly distinguish between these two alternative 
explanations. 

Pairbreaking effects are viewed as providing important in- 
sight into the origin of the cuprate pseudogap. Indeed, the dif- 
ferent pairbreaking sensitivities of T* and T c are usually pro- 
posed to support the notion that the pseudogap has nothing to 
do with superconductivity. To counter this incorrect inference, 
a detailed set of studies was conducted, (based on the BEC- 
BCS scen ario), of pairbreaking in the presence of impurities 
ll l4ll 1^1 and of magnetic fields IllfJl , These studies make 
it clear that the superconducting coherence temperature T c is 
far more sensitive to pairbreaking than is the pseudogap onset 
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FIG. 24: Comparison of theoretically calculated low T slope 
dHd/dT (main figure) for various doping concentrations (corre- 
sponding to different T c ) in the underdoped regime with experimen- 
tal data (inset) from Ref. llOSl The theoretical slopes are estimated 
using the low temperature data points accessed experimentally. The 
quantitative agreement is very good. 



temperature T*. Indeed, the phase diagram of Fig. 12 II which 
mirrors its experimental counterpart, shows the very different, 
even competing nature of T* and T c , despite the fact that both 
arise from the same pairing correlations. 



D. Anomalous Normal State Transport: Nernst Coefficient 

Much attention is given to the anomalous behavior of the 
Nernst coefficient in the cuprates ll45Tl . This coefficient is 
rather simply related to the transverse thermoelectric coeffi- 
cient a xy which is plotted in Fig. [25] In large part, the ori- 
gin of the excitement in the literature stems from the fact that 
the Nernst coefficient behaves smoothly through the supercon- 
ducting transition. Below T c it is understood to be associated 
with superconducting vortices. Above T c if the system were 
a Fermi liquid, there are arguments to prove that the Nernst 
coefficient should be essentially zero. Hence the observation 
of a non-negligible Nernst contribution has led to the picture 
of fluctuating "normal state vortices". 

The formalism of Ref.|47|can be used to address these data 
within the framework of BCS-BEC crossover. The results are 
plotted in Fig.[25]with a subset of the data plotted in the upper 
right inset. It can be seen that the agreement is reasonable. In 
this way a "pre-formed pair" picture appears to be a viable al- 
ternative to "normal state vortices". It will, ultimately, be nec- 
essary to take these transport calculations below T c . This is a 
project for future research and in this context it will be impor- 
tant to establish in this picture how superconducting state vor- 
tices are affected by the noncondensed pairs and conversely. 



V. CONCLUSIONS 

In this Review we have summarized a large body of work 
on the subject of the BCS-BEC crossover scenario. In this 
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FIG. 25: Calculated transverse thermoelectric response, which ap- 
pears in the Nernst coefficient, as a function of temperature for the 
underdoped cuprates. 



context, we explored the intersection of two very different 
fields: high T c superconductivity and cold atom superfluid- 
ity. Theories of cuprate superconductivity can be crudely 
classified as focusing on "Mott physics" which reflects the 
anomalously small zero temperature superfluid density and 
"crossover physics", which reflects the anomalously short co- 
herence length. Both schools are currently very interested 
in explaining the origin of the mysterious pseudogap phase. 
In this Review we have presented a case for its origin in 
crossover physics. The pseudogap in the normal state can be 
associated with meta-stable pairs of fermions; a (pseudogap) 
energy must be supplied to break these pairs apart into their 
separate components. The pseudogap also persists below T c 
in the sense that there are noncondensed fermion pair excita- 
tions of the condensate. 

The recent discovery of superfluidity in cold fermion gases 
opens the door to a set of fascinating problems in condensed 
matter physics. Unlike the bosonic system, there is no coun- 
terpart of Gross-Pitaevskii theory. A new theory which goes 
beyond BCS and encompasses BEC in some form or another 
will have to be developed in concert with experiment. As of 
this writing, there are four experiments where the simple mean 
field theory discussed in this review is in reasonable agree- 
ment with the data. These include the collective mode studies 
over the entire range of accessible magnetic fields |96, 97J]. In 
addition in the unitary regime, RF spectroscopy-based pair- 
ing gap studies l72l l9lll. as well as density profile l89ll and 
thermodynamic studies 17 3l 18811 all appear to be compatible 
with this theory. Interestingly, all of these provide indications 
for a pseudogap either directly through the observation of the 
normal state energy scales, T* and A, or indirectly, through 
the observation of noncondensed pairs. The material in this 
Review is viewed as the first of many steps in a long pro- 
cess. It is intended to provide continuity from one community 
(which has addressed the BCS-BEC crossover scenario, since 
the early 1990's) to another. 
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